Accepted for publication in The Astrophysical Journal Letters 2010 April 29 

Preprint typeset using I^-T^X style cmulatcapj v. 08/13/06 



O 

(N 



o 

u 

Oh- 
i 

o 

o3 



(N 
> 
O 

in 

O 
O 



STRONG GRAVITATIONAL LENSING BY THE SUPER-MASSIVE CD GALAXY IN ABELL 3827 

E. R. Carrasco 1 , P. L. Gomez 1 , T. Verdugo 2 , H. Lee 1 , R. Diaz 1 , M. Bergmann 1 , J. E. H. Turner 1 , B. W. Miller 1 

and M. J. West 1,3 

Accepted for publication in The Astrophysical Journal Letters 2010 April 29 

ABSTRACT 

Wc have discovered strong gravitational lensing features in the core of the nearby cluster Abcll 3827 
by analyzing Gemini South GMOS images. The most prominent strong lensing feature is a highly- 
magnified, ring-shaped configuration of four images around the central cD galaxy. GMOS spectro- 
scopic analysis puts this source at z ~ 0.2. Located ~ 20" away from the central galaxy is a secondary 
tangential arc feature which has been identified as a background galaxy with z ~ 0.4. We have mod- 
eled the gravitational potential of the cluster core, taking into account the mass from the cluster, the 
brightest cluster galaxy (BCG) and other galaxies. We derive a total mass of (2.7 ± 0.4) x 10 13 A/ Q 
within 37 h _1 kpc. This mass is an order of magnitude larger than that derived from X-ray observa- 
tions. The total mass derived from lensing data suggests that the BCG in this cluster is perhaps the 
most massive galaxy in the nearby universe. 

Subject headings: galaxies: clusters: individual (Abell 3827) — galaxies: elliptical and lenticular, cD 
— galaxies: evolution — galaxies: formation — gravitational lensing: strong 



1. INTRODUCTION 

As the densest galaxy environments known, the cores 
of massive clusters are expected to host the strongest 
dynamical evolution. Such cluster cores are found to 
be dominated by early-type, D or cD galaxies that are 
also the brightest cluster g alaxies (BCGs) . These have 
extended luminous halo es {S chombcrt 19 83), which are 
not necessarily smooth (j Johnstone et al .11 19911) and con - 



tain multiple or complex nuclei ( Rood fe Leirl I1979I ). 
They are located close to th e peak of the X-ray emis- 
sion (jJones fc Formanl Il984f ) and near the kinemati- 
cal centers of their clusters (jQuintana fe Lawriel 119821 : 
lOuintana et al.ll2000D . 

There has been a long-running debate over the 
extent to which BCGs can be assembled through 
continuous merging of galaxies in the cluster po- 
tential ( "galactic cannibalism" ) versus ea r ly merg- 
ing d uring cluster collapse (e.g. IWestJ 119941 : iDubinskil 
1998). However, recent work has shown that 
BCGs are likely to form via "dry" or dissipation- 



van DokkuirJ 120051 iBell et all 12006; 
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). These dry merg- 



ers are thought to be the primary mechanism through 
which massive galaxies grow from z ~ 1 to the present 
day and continue populating the upper end of the mass 
function without changing the overall mass density of el- 
liptical systems. It has been pointed out that the merger 
rate incre ases both with t he stellar mass of galaxies and 
with age (|Liu et al.ll2009t ). Therefore, there is a large 
fraction of major dry mergers in the nearby universe. In- 
deed, as ma ny as 3.5% of B CGs show ongoing evidence 
of mergers (|Liu et ai J 120091) . and about 10% - 20% of 
massive galaxies have undergone a dry merger in the last 
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gigayear dKhochfar fc Silkl [20091) . Thus, dry mergers ul- 
timately lead to the formation of very massive and dense 
BCGs with large mass-to-light ratios. 

We report here the discovery of a multi-component 
BCG that is, to our knowledge, the most massive galaxy 
ever seen in the local universe. This is located in the core 
of Abell 3827, a massive galaxy clusters in Abell's clus- 
ter catalog (richness class 2 and Bautz- Morgan type D, 
with an X-ray emission of L x (0.1-2.4 keV) = 2.1 x 10 44 
erg s _1 . This BCG is perhaps the most extreme exam- 
ple of ongoing galaxy cannibalism known: a super-giant 
elliptical that appears to be in the throes of devouring 
at least four other galaxies. Evidence for a recent ma- 
jor merger is also supported by the appearance of an 
extended asymmetric halo at the center of the cluster. 
The super-giant galaxy also shows features arising from 
strong gravitational lensing, the most prominent being a 
surrounding, magnified ring-shape configuration of four 
similarly-shaped images. The existence of such features 
provides a unique opportunity to study the mass distri- 
bution and the evolution of BCGs with unprecedented 
spatial detail. 

This Letter is organized as follows. In section 2, we 
summarize our observations and data reduction. In sec- 
tion 3, we describe our results, including the strong lens- 
ing modeling. In section 4 we present further discussion. 
Throughout this Letter we use a standard cosmology of 
H = 70 h km s" 1 Mpc" 1 , O m = 0.3 and f2 A = 0.7. 
At the redshift of Abell 3827, 1" corresponds to 1.83 h' 1 
kpc. 

2. OBSERVATIONS 

All imaging and spectroscopic data were collected 
with the Gemini Multi-Object Spectrograph (GMOS, 
iHook et all I2004T ) at the Gemini South telescope, in 
queue mode. 

The cluster was imaged with the <?', r' and i' filters 
on 2007 Nov 7 (UT), in photometric conditions, with 
good seeing (0'.'4 — 0'.'6). Observations were processed 
with the Gemini IRAF package (version 1.8) in the stan- 
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Fig. 1. — Color-composite image of Abell 3827 . The field of view is 5'x5' (~ 0.55 X 0.55 h -2 Mpc 2 at the distance of the cluster). The cD 
galaxy and its asymmetric halo are clearly seen at the center of the image. Top - left panel: the field - of - view is approximately 40"x40". 
The five nuclei contained in the the center are labeled N.l to N.5 and two stars are labeled S. Bottom-right panel: The field-of-view and 
orientation are similar to that of the top-left panel. The "central" arcs are labeled A.l to A. 4 and the tangential arc is labeled B.l. With the 
central nuclei removed and the image display scale chosen to maximize contrast, all of the gravitational-lensing features are more distinct. 



dard manner. Calibration on the standard magnitude 
system was ac hieved using observations of stars from 
lLandoltl (|1992). Object detection and photometry were 
performed on the r'-ban d image (2007 Nov 7 UT ) with 
the program SExtractor (|Bertin and Arnoutslll996ft . The 
MAG_AUTO parameter was adopted as the total mag- 
nitude for the detected objects. Colors were derived by 
measuring fluxes inside a fixed circular aperture of l'/5 
(~ 10 pixels) in all filters, corresponding to a physical 
aperture of 2.75 h _1 kpc at the distance of the cluster. 
All objects with SExtractor stellarity flag < 0.8 were se- 
lected as galaxies. We estimate that the catalog is com- 
plete down to r'=23 mag, since the number counts start 
to turn over at this value. The final catalog contains the 
total magnitudes, colors and structural parameters for 
747 galaxies brighter than r'= 23 mag. Figure Q] shows 



the color-composite g\ r' , i' image of the central region 
of Abell 3827 (0.55 x 0.55 /i~ 2 Mpc 2 ). Five central el- 
liptical galaxies (N.l to N.5) are embedded in a common 
asymmetric halo (top- left panel). The asymmetric halo 
has an ellipsoidal form, with a ~ 1.'3 and b ~ 1' for the 
semi-major and semi-minor axes, respectively. Strong 
lensing features detected around the central cD galaxy 
(at radii ~ 8" and ~ 20") are visible in the images lo- 
cated in the top- left and bottom-right panels of Fig. 1. 

Long-slit spectra of the central five elliptical galaxies 
were obtained on 2007 Jun 3 (UT) during bright time, 
through thin-cloud cover and ~ 1" seeing. The obser- 
vations were performed using the nod-and-shuffic mode, 
with the R400 grating centered at 8100A; this setup was 
chosen to minimize the effect of moon illumination and 
to obtain good sky subtraction. To avoid second-order 
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TABLE 1 

Main data of the five central ellipticals galaxies located in the core of Abell 3827 



Galaxy 


ID 


R. A. (2000) 


Decl.(2000) 


v M 


SV Z 


X 


y 


r' 


i' 








m s ^ 


(° ' ") 


(km s" 1 ) 


(km s -1 ) 


(arcsec) 


(arcsec) 


(mag) 


(mag) 


(mag) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


1323 


N.l 


22 01 53.99 


-59 56 45.14 


29548 


-115 


13.88 


-1.92 


17.03 


16.51 


0.51 


1324 


N.2 


22 01 53.31 


-59 56 43.14 


29763 


100 


3.75 


0.08 


16.19 


15.67 


0.52 


1325 


N.3 


22 01 52.69 


-59 56 41.22 


29951 


288 


-5.75 


2.00 


16.18 


15.56 


0.55 


1326 


N.4 


22 01 52.76 


-59 56 46.12 


28590 


-1073 


-4.63 


-2.91 


16.23 


15.81 


0.43 


1327 


N.5 


22 01 52.79 


-59 56 37.73 


33817 


4154 


-4.13 


5.48 


18.08 


17.70 


0.35 



Note. — (1) — Galaxy number in the Source Extractor catalog, (2) — Galaxy ID in Fig.l, (3) and (4) — R.A. and deck 
(J2000.0), (5) — heliocentric radial velocity. (6) — difference between the heliocentric radial velocities of each component 
and the systemic velocity, (7) and (8) — coordinates relative to the center of mass of the system located at ct2000 = 
22 h 01 m 53?06, fcooo = -59° 43' 43'.'22. (9) and (10) - total magnitudes, (11) - galaxy colors measured inside a fixed 
aperture of diameter 1'.'5. 



contamination from blue light in the red part of these 
spectra, the OG515 blocking filter was used. In addi- 
tion, multi-object spectroscopic (MOS) observations of 
galaxies in the Abell 3827 field were obtained on 2007 
Sep 7 (UT), during dark time, in photometric conditions 
and with seeing between 0'.'7 and 0'.'9. Two masks were 
observed with 0'.'75 slit widths, using the B600 grating 
centered at 5220A. During the observations, all spec- 
tra were dithered to account for the gaps in the GMOS 
CCDs. Spectra were reduced and calibrated in the stan- 
dard manner for each mode. 

3. RESULTS 

We have measured radial velocities for 67 galaxies us- 
ing cross-correlation techniques or emission-line fitting. 
Of these 67 galaxies, 56 are within ± 2500 km s _1 of the 
cluster redshift . Using the bi- we ight estimator for loca- 
tion and scale (jBeers at al.lll996T ). we calculate an aver- 
age velocity for the cluster of 29663 ±155 km s" 1 and 
a line-of-sight velocity dispersion of 1142 ± 125 km s _1 , 
with 55 member galaxies. The virial mass inside a ra- 
dius of 0.3 h- 1 Mpc is (3.861^1) x lO 14 /^ 1 M Q , with 
uncertainties at the 68% confidence intervals. 

In the cluster core, we confirm membership for four of 
the central elliptical galaxies. There are also two stars in 
the foreground (see Fig. 1). Three of the core galaxies 
have similar colors (see Table 1) and their velocities are 
consistent with the average cluster velocity (within <~ 300 
km s _1 ). Galaxy N.4 is slightly bluer than N.l, N.2 and 
N.3 and has a measured radial velocity ~ 1000 km s _1 
higher than the systemic velocity. Galaxy N.5 may not 
be in the cluster, but the background. 

3.1. Strongly lens ed features in the cluster core 

Thanks to the superb image quality delivered by 
GMOS and the Gemini South telescope, we have been 
able to detect several strongly-lensed features around the 
central cD galaxy in the core of the cluster. A thin, 
extended tangential arc (B.l in Fig. 1) is seen at ~ 
20" south-east from the central galaxy, with apparent 
length ~ 11" and thickness ~ 0'.'7. Surrounding the 
central core galaxies is a second and more prominent 
arc feature — a highly-magnified, ring-shaped configura- 
tion of four images around the central cD galaxy, with 
a radius of ~ 8". Only three of the arc components are 
visible in the inset color image in Fig. 1 (A.l, A. 2 and 
A. 3). The fourth image (A. 4 in Fig. 1) is between galax- 



ies N.l and N.2 and is obscured by the halo of the cD 
galaxy. All four features can be seen in detail in the 
bottom- right panel in Fig. 1 (see the caption for details). 
These features are ver y reminiscent of the lensed images 
seen in C10024+1654 (jCollev et al.lll996[ ). although the 
source galaxy in the present case is closer than z ~ 1 (see 
below) . 

The redshifts of the lensed sources were determined 
from long-slit observations taken on 2007 Nov 15-16, us- 
ing the R150 grating centered at 7150A. For the feature 
with the ring-shaped configuration (system A), we de- 
rived a redshift for one of the images located north-west 
of the cluster core (A. 2). Two emission lines {Ha and 
[N II]A6583A) are seen in the spectrum for A. 2; unfortu- 
nately no other lines were detected, due to light contam- 
ination from the cD halo. From the two emission lines, 
the estimated redshift for A.2 is z = 0.20443 ± 0.00073. 

The symmetric distribution of images in system A and 
the ring-shaped structure connecting arcs A.l, A.2 and 
A. 3 (and possibly A. 4) suggest that the images are of the 
same source galaxy. Therefore, we assume that all four 
images come from the same source and have the same 
redshift. This assumption should be taken with caution, 
without spectroscopy of the remaining system A images 
for confirmation, but is supported by our modeling, dis- 
cussed in the next section. 

The subtraction of light from the halo of the cD galaxy, 
as well as subtraction of all galaxies and the stars from 
the cluster core, highlights the ring-shaped structure con- 
necting arcs A.l, A.2, A. 3 and possibly A. 4 (see Fig. 1). 
The existence of this structure only ~ 15 h~ x kpc from 
the center of the cD galaxy gives us a unique opportunity 
to study the dynamics and mass distribution of a large 
BCG at unprcccdcntcdly small scales. 

Using several emission- lines ([Oil] A3727A, H(3, [OIII] 
AA4959,5007A) seen in the spectrum of the south- 
east tangential arc, B.l, we derived a redshift of z = 
0.40825 ± 0.00072. The relatively long, thin appear- 
ance of this arc suggests that the source galaxy lies very 
close to the inner fold caustic, assuming the lens clus- 
ter has an ellip tical spatial distribution of ma ss (e.g . see 
iFort k Meliien 11994 iNaravan k Bartelmannl Il999f ). If 
true, there must be three fainter counter-images: one 
to the north-east, one to the south-west and one at the 
center of the cluster core. 

3.2. Strong Lensing modeling and mass determination 
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TABLE 2 
Best-fitting parameters 



Comp. 


X 


Y 


e e 


Tcore 


fcut 


CO 




(") 


(") 


(°) 


(kpc) 


(kpc) 


(km s _1 ) 


Cluster+cD 


-0.9 


-2.5 


0.14 56 


32 


[1500] 


1174 


N.3(L*) 


[2.77] 


[2.82] 


[0] [0] 


[0.15] 


8 


256 



Note. Coordinates ar 

(a 2 ooo =22 h 01 m 53!15, feooo 
brackets are not optimized. 



relative to the center of image 2 
= -59° 56' 43'.'31). Values in square 



(2.7 ±0.4) x 10 13 Mq. This mass is slightly greater than 
that encl osed in Abell 1689 w ithin a similar radius (see 
Fig. 6 in iLimousin et al.ll2007|) . We estimated the mass 
of the cD galaxy as being ~ 30 % — 50% of the tota l 
mass of the cluster within 20" (jLimousin et al.1 [20071 ). 
This implies a cD galaxy mass between 8.1 x 10 12 M Q 
and 1.3 x 10 13 M©. Even assuming a conservative value 
of 30%, this would mean that the central cD galaxy in 
Abell 3827 is perhaps the most massive galaxy observed 
in the local universe. 



To produce the features described in the previous sec- 
tion, the cluster core must be very massive. In order to 
reconstruct the mass distribution in the core (the region 
enclosed by system A), we used the parametric model 
implemented i n the public l y ava ilable LENSTOOL 4 ray- 
tracing code (| Julio et al.l 120071) . To model the lens, 
we used two clumps of different scales — a large-scale 
halo, representing both the matter inside the giant cD 
galaxy and the dark matter halo of the cluster, and four 
smaller-scale clumps, modeling the small-scale perturba- 
tions associated with the galaxies embedded in the cD 
galaxy. For all components we adopt a dual Pseudo- 
Isothermal Elliptical Mass Distribution (dPIEMD; see 
lEhasdottir et al.ll2007t ). The dPIEMD can be character- 
ized by seven parameters: the center position (X, Y), the 
cllipticity e, the position angle 6 and the parameters of 
the density profile: the velocity dispersion <tq and two 
characteristic radii r core and r cut . For the large halo, 
we left most parameters free, with broad uniform priors. 
Only r cu t remained fixed at a value of 1500 h^ 1 kpc. 
In the small-scale clumps, the parameters of the density 
profiles were scaled as a funct ion of their galaxy lumi- 
nosities (see I Julio et all 120071 and references therm). 
Using as a scaling factor the luminosity L* , associated 
with the (/-magnitude of central galaxy N.3 (see Fig. 1 
and Table 1), we searched for the values of CTq an d Kut 
that yield the best fit, fixing r* ore at 0.15h _1 kpc. The 
best-fitting parameters of the model are listed in Table 
2. 

Because the ring structure shows multiple subcompo- 
nents, we selected the most reliable set of multiple im- 
ages and associated them with three different background 
sources in order to make our fit. Considering an un- 
certainty in the position of any image equal to 0'.'2, we 
found after the optimization in the image plane a Xdof 
= x 2 /DOF = 69/4 ~ 17. Our model reproduces well the 
positions of the observed subcomponents, with a mean 
scatter less than 0'.'8. To highlight this, we show in Fig- 
ure [5J the model-predicted counter-images of the sub- 
structures in arc A.l (green crosses), which are in agree- 
ment with the image positions for all arcs (orange cir- 
cles). The fit also predicts a central demagnified image 
which is lost in the cD light distribution. It is impor- 
tant to note that our model is oversimplified, since we 
are assuming spherical halos for the galaxies. A detailed 
model of Abell 3827 is beyond the scope of the present 
work. 

Given our best model (see Table 2), we calculated 
the total mass inside a radius of 20" (location of the 
tangential arc B.l at ~ 37 h _1 kpc) and found M = 

4 http:/ /www. oamp.fr/cosmology/lenstool/ 
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Fig. 2. — g'- band image of Abell 3827 with the local median 
average subtracted. North is up and east is left. The external 
critical line (red) and the associated caustic line (gray) for a source 
at z = 0.2 are depicted. Orange circles are the measured image 
positions. Green crosses represent the predicted image positions, 
using as input the three substructures that constitute arc A.l. We 
also show with small asterisks (cyan) the positions in the source 
plane associated with these three points. For clarity, we do not 
depict the demagnified images produced by strong galaxy - galaxy 
lensing events. 

4. DISCUSSION 

Using the strong-lensing model presented above, we 
have analyzed the mass distribution to unprecedented 
spatial resolution. The cluster core is spatially very con- 
centrated. Even if we assume a conservative fraction for 
the mass of the central cD galaxy (see above) , this galaxy 
is very massive. The cD galaxy could be an extreme 
example of the effects of dry mergers on the mass of 
BCGs; dry mergers can produce an increase (by a factor 
of up to 3) in the dark-matter-to- stellar mass ratio for the 
most massive systems at present (jRuszkowski fc Springell 
[2001 . 

We have compared the total mass derived from our 
strong-lensing analysis with that derived from the X-ray 
gas. We have calculated the total mass as a function 
of radius for Abell 3827, based on a recent X-ray anal- 
ysis by L. Valkonen et al. (2010, private communica- 
tion) . We derive the total cluster mass by assuming 
hydrostatic equilibrium and using their best fit for the 
three-dimensional distribution of X-ray emitting gas. L. 
Valkonen et al. (2010, private communication) report 
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radius (Mpc) 

Fig. 3. — Mass profile derived from the best fit model for the X- 
ray gas distribution (continuous line). The dash-dotted line shows 
the profile of the X-ray distribution for a more typical cluster. 
The short-dashed line shows the mass profile derived from strong 
lcnsing. The vertical solid line shows the location of the tangential 
arc (B.l) at 37 h" 1 kpc. 

M 2 oo = 9.7 x 10 14 M inside a radius r 20 n = 2.84 h" 1 
Mpc, derived from the M — T relation bv lSalhen et al.l 
(|2009t) for a cluster temperature of kT = (7.15 ± 0.18) 
keV. Fig. [3] shows the radial distribution of the total 
mass derived from the X-ray gas and from the strong- 
lcnsing analysis. For the X-ray gas, we plot mass pro- 
files: (1) for a typical cluster with r core — 250 h _1 kpc 
and j3 = 0.7 (dash-dot line), and (2) for the derived 
total mass of Abcll 3827 within 37 h" 1 kpc (solid ver- 
tical line). Finally, we plot the mass profile for Abell 
3827 derived from strong lensing as a dashed line. At 
all radii, the mass derived from strong lensing is at least 
a factor of 10 larger than that derived from X-ray data. 
Discrepancies between strong-lensing cluster masses and 
X-ray cluster masses for the same clust ers have been re- 
ported previously in the literature (e.g. lGitti et al.ll2~007l : 
lHalkola et al.l l2008h . Moreover, significant biases of up 
to ~ 50% may be introduced into strong-lensing mass 
estimates when models are extrapo lated outside the Ein- 
stein ring (|Meneghetti et all 12003). Analysis using both 
strong-lensing and X-ray methods at small radii show 
that dif ferences can grow by as much as an order of mag- 
nitude (jVerdugo et al.ll2007l) . 

How can we explain this discrepancy? For the lower 
limit on mass, within 3cr error of the strong-lensing fit, 



we obtain 6.3 x 10 12 M© inside the ring structure of ra- 
dius 10" (~ 19 h _1 kpc). Hence, the most likely un- 
certainties, related to extrapolation of the strong-lensing 
model to larger radii, are unable to explain our factor- 
of-ten discrepancy. Assuming a spherical mass distribu- 
tion, the mass inside system A is (9.6 ± 0.9) x 10 12 M©. 
The large difference between the strong-lensing mass and 
X-ray mass could be related to the dynamical state of 
the cluster, to the mass inferred from the X-ray gas, 
or to both effects. Total X-ray cluster masses derived 
from flat-core density profiles may be underestimated by 
at least a factor of 2 within the cent ral ~ 30 h -1 kpc 
(IVoigt fc Fabian! 120061) . Simulations (|Meneghetti et al.l 
120091 ) of strong-lensing models have shown that physi- 
cal substructures along the line of sight to the cluster 
can overestimate the derived total mass by a factor of 2; 
in fact, our GMOS spectra reveal a bi-modality in the 
velocity distribution of cluster galaxies, suggesting that 
Abell 3827 is presently merging. In all, these effects can 
produce a factor of about 4, still insufficient to explain 
the large discrepancy between strong-lensing and X-ray 
derived masses. Moreover, if we take out the factor of 4 
due to the two effects described above, we are left with a 
cD galaxy with a mass of ~ 2 x 10 12 M©, which appears 
to be one of the most massive galaxies known in the local 
universe. Detailed studies of the central objects at differ- 
ent wavelengths will reveal insights into the origins and 
subsequent evolution of this peculiarly massive galaxy. 
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